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Abstract: The lowest radiative correction to the Casimir energy density between two parallel 
plates is calculated using effective field theory. Since the correlators of the electromagnetic field 
diverge near the plates, the regularized energy density is also divergent. However, the regularized 
integral of the energy density is finite and varies with the plate separation L as 1/L 7 . This 
apparently paradoxical situation is analyzed in an equivalent, but more transparent theory of 
a massless scalar field in 1+1 dimensions confined to a line element of length L and satisfying 
Dirichlet boundary conditions. 
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£T) ' The remarkable prediction of quantum electrodynamics (QED) by Casimir|IJ that 

two parallel metallic plates in vacuum and with separation L should be attracted by a 
force per unit area equal to 7r 2 /240L 4 , has very recently been verified experimentally 
by Lamoreaux^] with much higher precision than before||. In fact, the experiment is so 
accurate that with further improvements it can hopefully also demonstrate the corrections 
due to finite-temperature effects j§] |5| fj] [0 . It would then be a new macroscopic system 
Sh I where the intricacies of quantum field theory can be studied. 

q-( Despite an enormous literature on the subject (see for example the review article^] 

and a recent textbook^]) very little has been done to investigate the radiative corrections 
to the this effect based on QED which is basically an interacting theory. It is still not clear 
at which order in the fine structure constant a the first correction will appear and how 
it will vary with the plate separation. One of the first calculations were done by Bordag, 
Robaschik and Wieczorek[10] who found a correction proportional to a/L 5 in lowest order 



QED perturbation theory using a photon propagator which is appropriately modified due 
to the plate confinement. Using instead periodic boundary conditions, Xue found that the 
corresponding radiative corrections were exponentially small [ffl 



Instead of using QED, we will here use the modern approach of effective field theory 
JT^ ]. (A pedagogic introduction has been given by Kaplan fl3j]). The energy scale relevant 
for the Casimir effect is much smaller than the electron mass m. We can thus integrate 
out the heavy electron modes from the full Lagrangian and thus obtain a low-energy 
effective Lagrangian for the interacting photon field. The interactions are due to coupling 
of the photon to virtual electron-positron pairs in the vacuum and they become local in 
the low-energy limit. To lowest order in a the first effective coupling is just the Uehling 
correction^] which modifies the photon propagator at short distances. It gives rise to 
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the effective photon Lagrangian 

where F^ u = d^A u — dyA^ is the electromagnetic field strength. If this new interaction 
contributes, we see that it will modify the Casimir energy density by a term proportional 
with a/m 2 L 6 . Barring some extra factors of m from the evaluation of the corresponding 
Feynman diagram or from the integration over the volume between the plates, we see that 
we are not able to reproduce the result of Bordag et al [p^Of] . However, the corresponding 
Feynman diagram which is a photon loop with the Uehling correction in the propagator, is 
in fact zero after regularization. This is most easily seen by using the equation of motion 
OFuv = for the photon following from the leading, Maxwell part of the Lagrangian 
(II). One might perhaps argue that the classical equation of motion only applies to free, 
i.e. on-shell photons and not to virtual particles inside a Feynman diagram. But a more 
careful analysis |l5| shows that the Uehling interaction does not contribute to the vacuum 
energy. 

The lowest order correction must thus arise from higher order interactions, i.e. of 
dimensions eight or more. Among these we first have the next order Uehling term propor- 
tional with F^U 2 . But this will again not contribute for the same reasons as above. 



But to the same order we also have the Euler-Heisenberg interaction[16] 
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(E^ — B 2 ) 2 + 7(E ■ B)^ (2) 



which describes the interaction of four photons via coupling to a virtual electron loop 
in the vacuum. It has previously been used by Scharnhorst and Barton to investigate 



the propagation of light in the same geometry of two parallel plates [17]. For dimensional 
reasons it will thus give an additional Casimir force which varies with the plate separation 
like 1/L 8 . 

In lowest order perturbation theory the correction to the vacuum energy density is 
given by the expectation value ( C eh ) ■ For its evaluation we need the correlators ( E 2 ) 
and ( B 2 ) which give the fluctuations of the free electric and magnetic fields between the 
plates. The modes must satisfy the metallic boundary conditions nAE = n-B = at the 
plates and have previously been obtained by Liitken and Ravndal by quantizing the field 
in the Coulomb gauge and expanding it in electromagnetic multipoles|18|]. They can also 
be obtained from the coincidence limit of the photon propagator constructed by Bordag 
et al |l0| in the Lorentz gauge JIB]]. The results can be summarized as 



where = zir/L is the scaled distance between the plates and the function 

Id 3 ,„ 3 2 
sin tt sin 



F ^ = -2W COte = —»-—- (4) 
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gives the position dependence of the fluctuations. However, when one calculates the energy 
density £ = ^(( E 2 ) + ( B 2 }), it cancels out and gives £ = — 7r 2 /720L 4 which is the standard 
result for the Casimir energy per unit volume Q . When one gets very near one of the plates, 
the fluctuations become large. Near the left plate at z = they increase like 

This represents also the fluctuations near any curved surface since it will always look plane 
when one gets near enough [p!s|] . 

With the free field electromagnetic correlators established between the plates, it is 
now straightforward to calculate the lowest order correction from the Euler-Heisenberg 
interaction (^). In terms of Feynman diagrams, it is simply given by the contribution 
from the product of two photon loops, 
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This additional term in the Casimir energy density is now seen to diverge when one 
gets near one of the plates. The divergence is directly related to the the corresponding 
divergence in the correlators. We are now faced with the problem that when calculating 
the resulting correction to the Casimir force, we need the total vacuum energy between 
the plates which is obviously also divergent. 

This is physically meaningless, but not really a new situation. For instance, when 
one calculates the regularized vacuum energy density for free photons inside a spherical 
cavity [^], it diverges near the surface as in (||) and the integrated energy is thus also 
divergent. However, if one instead calculates the total energy which afterwards is regu- 



larized, one obtains a finite result [21] [22] [23]. That the processes of regularization and 



integration do not commute in general for these kinds of problems have been discussed by 



Deutsch and Candelas[19| 



In our case with interacting fields we can obtain the total unregularized energy by 
integrating the unregularized energy density. The resulting expression can then be reg- 
ularized and we have a finite, total Casimir energy. Since the calculation is somewhat 
cumbersome, we can illustrate it by neglecting the uninteresting degrees of freedom par- 
allel to the plates and consider instead the simpler theory of a massless field in 1+1 
dimensions with Lagrangian £ = (1/2) [(c^0) 2 — (<9 2 0) 2 ]. Here dt4> can be taken to be 
the electric field and d z 4> to be the magnetic field. We can thus continue to talk about 
the electric £e = (l/2)( (dt 4>) 2 ) and magnetic £b = (l/2)( {dt4>) 2 ) contributions to the 
vacuum energy density. Imposing the Dirichlet condition that the field vanishes at the 
boundaries, the normalized eigenmodes of the field are 

T 

- sin (u n z) (7) 
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where oo n = irn/L with n = 1, 2, 3, ... . The total vacuum energy is therefore 

oo 

Eo = ^y n = -— (8) 

71=1 

when we use zeta-function regularization where ((— 1) = —1/12. Assuming a constant 
vacuum energy density, it is therefore u = — 7r/24L 2 . 

Let us now consider the partial contributions to the energy. The part coming from 
fluctuations in the electric field is given by the divergent sum 

^ oo oo 

£ E = -^Y^^n sin 2 (u n z) = — Y ^„(1 - cos 2u n z) (9) 

n=l n=l 

If we integrate over the 1-dimensional volume to obtain the total Casimir energy, the first 
term gives — 7r/48L after zeta-function regularization. But in the last position-dependent 
part, we see that each term gives zero because of the imposed boundary conditions. We 
are thus left with a finite contribution to the total energy which is exactly one half of the 
full energy. The other half comes from the magnetic contribution. 

Let us now instead derive a finite result for the energy density. It can be written as 
where 

oo 

S(0) = £sin20n (11) 

71=1 

is a divergent sum. Using again zeta-function regularization, it becomes S reg {6) = ^cot# 
and we have the regularized energy density 

e * 9 = ~m? (3 " ikTfl) (12) 

This corresponds to the result in (|3|) giving the fluctuations in the electric field between 
plates. The last term gives the dependence on the position in the volume. It is positive 
definite and diverges when we approach the boundaries of the volume, i.e. the end points 
of the interval. As a consequence the corresponding integrated energy is also infinite. 
This is in sharp contrast to the result above where the integrated contribution from the 
position-dependent term gave zero. 

The underlying problem is obviously the behaviour of the fluctuations near the bound- 
aries where sin 9 = 0. Then the function ([ll]) is zero while S reg is infinite. In order to better 
understand this apparent paradox, we can use a more physical regularization method based 
upon an exponential cutoff instead of the previous more mathematical approach using the 
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analytical continuation of the zeta-function. We thus define the sum (|ll| ) by the limit 
S(9) = lim e _+o S(e, 6) of the convergent sum 



S( e , 9) = Y e~ m sin 29n = 6 - ( 1 :', ) 

v ; ^ 1 - 2e~ e cos 29 + e~ 2e 

n=l 



Since it is the derivative of this sum which appears in the energy density @, it is seen 
not to contribute to the integrated energy even with a non-zero cutoff. This is consistent 
with what we found using zeta-function regularization. However, the regularized energy 
density is found in the limit e — > from 

S(e,9) = 1 -cot9- 1 -^-e* + 0(e*) (14) 
2 o sin 9 

As long as sin# > 0, i.e. inside the boundaries, we can take the e = and we recover 
the zeta-function result. But when we get sufficiently close to the boundaries, the ex- 
pansion breaks down and we have in general a cutoff-dependent result. This should not 
come as a complete surprise since the energy density on the boundaries must be strongly 
dependent on the microscopic and physical properties of the material in the boundaries. 
But the welcome result is that the integrated energy is independent of the cutoff and the 
regularization method. 

One can emulate the Euler-Heisenberg interaction in our one-dimensional system by 
considering the Lagrangian 



2 m z 



Here a is some small dimensionless constant and m a heavy mass. The Lagrangian is in- 
variant under the field transformation <p — » (f) + const and thus describes massless particles. 
Treating the interaction in lowest order perturbation theory, one easily finds the resulting 
vacuum energy density after regularization to be 



24L 2 8m 2 L 4 \18 sin 4 

Now it diverges near the endpoints where sin# — > and the integrated Casimir energy is 
infinite. However, if we instead first integrate the energy density and then regularize as 
above, we get the finite result 

2 

7T OiTT 

E ° = "24L " 144m 2 L 3 (17) 

The potentially divergent part now vanishes since £(— 2) = 0. We see that the total 
Casimir energy is furnished by just the constant part of energy density (|^). Again there 
is no contribution from the position-dependent terms. 

With this approach we can now calculate a finite, total Casimir energy for the photons 
between two plates and interacting via the Euler-Heisenberg interaction (pi). Zeta-function 
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regularization must now be combined with dimensional regularization for the transverse 
degrees of freedom. As a result we hnd again that the position-dependent part of the 
energy density @ does not contribute to the integrated energy which now becomes [|l5|] 

E = - (18) 

720L 3 2 7 3 5 5 3 m 4 L 7 1 ' 

Since the energy density between the plates can thus effectively be taken to be constant in 
calculating the total energy, it could also have been obtained more directly by requiring 
instead the field to be periodic between the plates with period (5 = 2L. This will then 
give a constant energy density £ = Eq/L also for the interacting field. Letting L — » 1/2T 
in the expression for £ we thus have the free energy density for an interacting photon 
gas in thermal equilibrium at temperature T. The result is in agreement with a previous 



calculation by Barton using a semi-classical method [25]. 
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